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We study nonlinear bound states, or solitary waves, in the Dirac-Maxwell system proving the existence of solu- 
00 ' tions in which the Dirac wave function is of the form (f>(x, o;)e _l " t , ui 6 (— m, w*), with some ui* > — m, such that 

(N \ 4>{-M G H 1 (R S , C 4 ), ,oj)\\ 2 L 2 = 0(m - |w|), and \\<j>( • , u)\\ L °° = 0(m - |w|). The method of proof is 

an implicit function theorem argument based on an identification of the nonrelativistic limit as the ground state of the 
Choquard equation. 

O-i 

*£h ■ 1 Introduction and results 

The Dirac equation, which appeared in [Dir28] just two years after the Schrodinger equation, is the correct Lorentz 
invariant equation to describe particles with nonzero spin when relativistic effects cannot be ignored. The Dirac 
equation predicts accurately the energy levels of an electron in the Hydrogen atom, yielding relativistic corrections 
to the spectrum of the Schrodinger equation. Further higher order corrections arise on account of interactions with 
the electromagnetic field, described mathematically by the Dirac-Maxwell Lagrangian, which aims to provide a self- 
consistent description of the dynamics of an electron interacting with its own electromagnetic field. The perturbative 
£N| ■ treatment of the Dirac-Maxwell system in the framework of second quantization allows computation of quantities 

such as the energy levels and scattering cross-sections, which have been compared successfully with experiment; of 
course this quantum formalism does not provide the type of tangible description of particles and dynamical processes 
familiar from classical physics. Mathematically, the quantum theory (QED) has not been constructed, and indeed 
may not exist in the generally understood analytical sense. In particular it is a curious fact that although the electron 
is the most stable elementary particle known to physicists today, there is no mathematically precise formulation and 
proof of its existence and stability. This has resulted in an enduring interest in the classical Dirac-Maxwell system, 
both in the physics and mathematics literature. Regarding the former, the relevance of the classical equations of 
motion for QED has been widely debated. The prevalent view seems to be that the Dirac fermionic field does not 
have a direct meaning or limit in classical physics, and hence that the classical system is not really directly relevant 
to the world of observation. Nevertheless, there have been numerous attempts, both by Dirac himself and by many 
others - see IIDir621 IWak66l |Lis951 and references therein - to construct localized solutions of the classical system or 
some modification thereof, with the aim of obtaining a more cogent mathematical description of the electron (or other 
fundamental particles). 

We consider the system of Dirac-Maxwell equations, where the electron, described by the standard "linear" Dirac 
equation, interacts with its own electromagnetic field which is in turn required to obey the Maxwell equations: 

[^(idp - eA^ij) -rm/j = 0, 
yd^dpA" = J v , 8^ = 0, 

with the charge-current density J M = (p, J) generated by the spinor field: 

,P l = eVn"V- (1-2) 

Above, p and J are the charge and current respectively. We denote t/j = ("f ip)* = "0*7°, with ip* the hermitian con- 
jugate of ip. The charge is denoted by e (so that for the electron e < 0); the fine structure constant is the dimensionless 
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coupling constant a = i- « 1/137. We choose the units so that h = c = 1. We have written the Maxwell equations 
using the Lorentz gauge condition d^A^ = 0. The Dirac 7-matrices satisfy the anticommutation relations 

{7 M ,7 v }=2<r, 

with g' 1!y = diag[l, — 1, —1, — 1]. The four- vector potential has components (A , A), with A = (A 1 , A 2 , A 3 ), so 
that the lower index version A^ = g^A 1 ' has components (A , —A) so Aq = A . Following 1BD64I and 1BS77I . we 
define the Dirac 7-matrices by 

= (A i) ■ v - (* _y . 

where I2 is the 2x2 unit matrix and a, are the Pauli matrices: u\ = ( ? _ | , cro = ( ? „ | , 03 = f _ ^ 

yi 0; V 1 °/ V° - 1 

After introduction of a space time splitting, the system ( II. U takes the form 

id t if) = Of(-iV-eA)il} + m^ + eA if), (df - A)A° = etp*ifi, (<9 2 - A)A = eip*aip. (1.4) 
Here a = (a 1 , a 2 , a 3 ), and a- 7 and /3 are the 4x4 Dirac matrices: 



(1.5) 



h 
-h 

with {(Jj}j =1 the Pauli matrices. We will not distinguish lower and upper indices j of a and a, so that ctj = a 3 , 
(jj = cr- 7 . The a-matrices and 7-matrices are related by 

7 J '=/3a J , l<i<3; 7 = /3. 

Numerical justification for the existence of solitary wave solutions to the Dirac -Maxwell system (II . lb was obtained 
in HLis95l . where it was suggested that such solutions are formed by the Coulomb repulsion from the negative part of 
the essential spectrum (the Klein paradox). The numerical results of QLis95l showed that the Dirac-Maxwell system 
has infinitely many families of solitary wave solutions (f>N{x,uj)e~' luJt , ui > —m. Here the nonnegative integer N 
denotes the number of nodes of the positronic component of the solution (number of zeros of the corresponding 
spherically symmetric solution to the Choquard equation; see iQ. A variational proof of existence of solitary waves 
for ui G (—m, 0) and with N = first appeared in |EGS96|, and the generalization to handle ui € {—ni, m) is in 
IIAbe98l . 

In the present paper, we give a proof of existence of solitary wave solutions to the Dirac-Maxwell system based 
on the perturbation from the nonrelativistic limit and also obtain the precise asymptotics for the solution in this limit. 
The physical significance of these types of solitary wave solutions requires not only their existence but also stability, 
and it is to be hoped that the type of detailed information about the solutions which is a consequence of the existence 
proof in this article, but does not seem to be so easily accessible from the original variational constructions, will 
be helpful in future stability analysis (see Remark [3~6l below) . The second motivation for presenting this proof is 
to realize mathematically the physical intuition explained in BLis95l which explains the existence of these bound 
state solutions in terms of the Klein paradox ( HBD64I §3.3]). Moreover, once one knows that the excited eigenstates 
of the Choquard equation are nondegenerate (currently this nondegeneracy is established only for the ground state, 
N = [Len09|), our argument will yield the existence of excited solitary wave solutions in Dirac-Maxwell system, 
extending the results of [EGS96| to N > 1. We will construct solitary wave solutions by deforming the solutions to 
the nonrelativistic limit (represented by the Choquard equation) via the implicit function theorem. Such a method was 
employed in lOjinOO Gua08 .| for the nonlinear Dirac equation and in IIRNlOal IStulOl IRNlObl for Einstein-Dirac and 
Einstein-Dirac-Maxwell systems. 

The solitary wave ((j>e~ tuJt , A ti (x)) satisfies the stationary system 

u><t> = ol ■ (-iV - eA)0 + m/3<t> + eA°(f>, -AA ft = e^(f>. (1.6) 
Theorem 1.1. There exists w* > —ni such that for to G (~m, w*) there is a solution to il.6i of the form 



e 3 <?i(ex,e) 
e 2 $ 2 (ex,e) 
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with 

G C°°((0,e*); (H 1 {R 3 ;C 2 )®H 2 (M 3 ;C 2 ))), e* = ^m 2 -^, 
and with 

A* e C°°((0,e*) ; iJ^R 3 ^) n L°°(R 3 ,R)), < /i < 3. 

Above, H 1 = i? 1 (IR 3 ,K) is the homogeneous Dirichlet space of L 6 functions with \\f\\ 2 ^ 1 = J |V/| 2 cfe < oo. For 
smaW e > 0, owe /za.s 

||*i - «Mh* = 0(e 2 ), ||* 2 - *a|| H i = 0(e 2 ), 

= lim#i(j/,e), * 2 (y) = lim* 2 (y,e) 
are of Schwartz class. The solutions could be chosen so that in the nonrelativistic limit e = one has 

*a(y) = Po(l/)n, #i(y) = ^-a-V y £ 2 {y), (1.7) 

where n G C 2 , |n| = 1, one/ <^o w a strictly positive spherically symmetric solution of Schwartz class to the Choquard 
equation 

~ 2^° = ~2^ A<A) ~ (idd *^")^ ' <?o(*)eR, iel 3 . (1.8) 

Remark 1.2. The existence of a positive spherically-symmetric solution c/5q G ,5^(R 3 ) to ( 11.81 ) was proved in IILie77l . 

Here is the plan of the paper. We give the heuristics in §2] The Choquard equation, which is the nonrelativistic 
limit of the Dirac-Maxwell system, is considered in <|3] In [JH we complete the proof of existence of solitary waves 
via the implicit function theorem. 

2 Heuristics on the nonrelativistic limit 

The small amplitudes waves constructed in Theorem ll.il are best understood physically in terms of the non-relativistic 
limit. Since we have set the speed of light and other physical constants equal to one, the relevant small parameter is 
the excitation energy (or frequency) as compared to the mass m. To develop some preliminary intuition regarding the 
non-relativistic limit, following IILis95L we neglect the magnetic field described by the vector-potential A 3 , getting 

id t ip = -ia-Vtp + m/3ip + eA°ip, {d 2 - A)A° = e^*i> . 

Let us consider a solitary wave solution t()(x,t) = 4>(x)e~ luJt , with <p(x) = 

A a = A°(x) only. Then 4>u 4>2, and A satisfy 

(u - m)0i = -ier-V0 2 + eA°<f>i, (to + uj)(j) 2 = -ia-Vfa + eA°<f) 2 , -A A = e($>i + <i>* 2 <i> 2 ) , (2.1) 

where a = (oi, cr 2 , 0-3), the vector formed from the Pauli matrices. Consider small amplitude solitary waves with 
uj ps —to. Then A Q is small and — 2to</>i ps — £<x-V<^>2> 

(to + uj)4> 2 ps —Acf) 2 + eA°<p 2 , -AA" ps e(0*0i + (f> 2 cj) 2 ). 

ZTO 

Denoting e 2 = m 2 — w 2 , < e <C to, the above suggests the following scaling: 

V = ex, d x = ed y , A°(x) = e 2 A°(ex), <t> x (x) = ^lN, M*) = ^ 2 (ex). (2.2) 

Note that since <f>j and A° depend on uj and x, the scaled functions A and <Pj are functions of y and of e. In the limit 
e^O, denoting 

<£=lim<Z>, A° = limA°, 

we arrive at 

- 2m$x = -ia-V y £ 2 , -^-£ 2 = -icr-VyS^ + eA°£ 2 , -A v A° = e£* 2 <P 2 , (2.3) 




b 2 (x) 



where <b\, 6 2 G C and 
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which can be rewritten as the following equation for <&2 only: 



1 -&2 = -—A y S 2 - eA°# 2 , 



-A y A° = e£*&' 



2^2 , 



(2.4) 



2m 2m 

with the understanding that <Pi is then obtained from the first equation of (12.3b . 

Remark 2.1. Regarding self-consistency of this approximation: one can check that, when using the scaling ( 12.2l i. the 
magnetic field vanishes to higher order in the limit e — s- 0, in agreement with 0Lis95l . Indeed, A = — A -1 J, where 
J = eip*aip = 0(e 5 ), hence A = — A -1 J = 0(e 3 ). The second equation from ( 12. Q would then take the form 

(m + u>)(f)2 = — icr-V^i — eA-cr(f)i + eA°02, 

where A-cr<fri = 0(e 6 ) while other terms are 0(e 4 ). Thus the approximation is al least self-consistent, and the 
analysis in SjUjustifies this rigorously. 

Remark 2.2. Regarding symmetry: while it is clear that radial symmetry of both <\>\ and (f>2 is inconsistent with (12 .3b - 
solutions of the form 



g( r ) 



cosy 
1 



if(r) 



cos 9 
-e"^ sin 9 



(2.5) 



are permitted in principle, suggesting that in the non-relativistic limit <p2 could be radial, or to be more precise of the 



form <p2{ll) = 







<G C 2 , where the spherically symmetric function $(y) € C is to satisfy 



2m 2m v 



-A y A» 



el*l 



(2.6) 



The starting point for our perturbative construction of solitary wave solutions to ( 11.41 ) is indeed a radial solution of 
( 12.6b , although the exact form of these solitary waves has to be modified from ( 12.5b when the effect of the magnetic 
field B = V x A is included, see |Lis951 §5]. The method of proof we employ does not require any particular 
symmetry class of the solitary wave. 

The above discussion suggests that the system (12.6b will determine the non-relativistic limit to highest order. 
The system d2.6b describes a Schrodinger wave function with an attractive self-interaction determined by the Poisson 
equation. Because the sign of the interaction is attractive it is often referred to as the stationary Newton-Schrodinger 
system. It is equivalent to a nonlocal equation for $ known as the Choquard equation, which is the subject of the next 
section. 



3 The nonrelativistic limit: the Choquard equation 

The system ( 12.6b can also be obtained by looking for solitary wave solutions in the system 
1 



id t i> 



2m 



-AV = eip*i/>, ip(x,t)eC, V(x, t) G M, x € 



(3.1) 



This is the time-dependent Newton-Schrodinger system. If [6e V(x)) is a solitary wave solution, then <f> and V 
satisfy the stationary system 



OJtp 



2m 



Acj) - eV<p, 



-AV = e|0| : 



We rewrite the system ( |3.1b in the non-local form, called the Choquard equation: 



idti> = + e^GVf )V, 1>(x,t) e C, 

2m 



x e 



(3.2) 



(3.3) 



where A 1 is the operator of convolution with — j^jjj ■ The solitary waves are solutions are of the form ?p(x, t) = 
<f) ul {x)e~' tu -' t , with W satisfying the non-local scalar equation 

1 



UJ(j) 



2m 



A^ + e^- 1 ^! 2 )^. 



(3.4) 
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This suggests the following variational formulation for the problem: find critical points of 

£ Choquard (0) = ^ / -£-ff lH ^ m f r dxdy , (3.5) 

2m J &TrJJ \x-y\ 

subject to the constraint J \cf>(x)\ 2 dx = const. This formulation is the basis of the existence and uniqueness proofs 
in the references which are summarized in the following theorem. 

Lemma 3.1 ( 1Lie77l ILio80l IMZ10I ). For all lo < and N e Z, N > 0, the equation (13.3b admits solitary wave 
solutions 

il>(x,t) = <pn(x, Lo)e~ lu)t , lim ipjy(x 7 uj) = 0, 

\x\— ^oo 

with lpm{x 1 lo) a spherically symmetric solution of (13. 4t . These solutions differ by the number N of zeros (or nodes), of 
the profile functions u?n{x, lo), considered as a function ofr = \x\. The profile function ipo with no zeros minimizes the 
value of the energy functional Echoquard(4>) amongst functions with fixed L? norm, and is the unique (up to translation) 
positive H 1 solution of d3.41 >; the corresponding solitary wave is known as the ground state. 

Remark 3.2. Together with the heuristics in the previous section, this result suggests that for u sufficiently close to 
— m there might exist infinitely many families of solitary waves to the Dirac-Maxwell system, which differ by the 
number of nodes. As mentioned in [EGS96], the variational methods used in that paper are hard to generalize to prove 
the existence of multiple solitary waves for each lo (such a multiplicity result is obtained in [EGS96 | for the Dirac - 
Klein-Gordon system). 

Remark 3.3. The <fr(x) and V(x) for different values of uo < can be scaled to produce a standard form as follows. 
Let C > satisfy £ 2 = — lo and write y = (x , <f>(x) = ( 2 u((x) , and V(x) = C, 2 v(C,x). Then (13. 2t is equivalent to the 
following system for u(y), v(y): 

— u = — - — Ayii — evu, — A y v = e\u\ 2 . (3.6) 

In the remainder of this section we summarize the properties of the linearized Choquard equation which follow 
from [Len09 | and are needed in 2] Consider a solution to the Choquard equation of the form 

ip{x, t) = {ipo(x) + R{x, t) + iS{x, t))e- lU0 \ 

with R, S real-valued. The linearized equation for R, S is: 



di 



R 







Lo 




~R 


S 











S 



Lo = -—A-LOo + e'-A-H^l), Li = L + 2e 2 A" 1 (^ • )<^o- (3.7) 
Zm 



Both Lq and L\ are unbounded operators L 2 — > L 2 which are self-adjoint with domain H 2 C L 2 . Clearly Lotpo = 0, 
with G ad(Lo) an eigenvalue corresponding to a positive eigenfunction ipo\ it follows that is a simple eigenvalue 
of Lq, with the rest of the spectrum separated from zero. The range of Lq is {ipo}' 1 , the L 2 orthogonal complement of 
the linear span of <po. Notice that L\ = ^(E"(tp ) — looQ" '(<fo)) ■ 

Lemma 3.4. The self-adjoint operator L\ : H 2 — > L 2 has exactly one negative eigenvalue, which we denote — Ao, 
and has a three dimensional kernel Ker L\ spanned by {dj(po}j—i- The range of L\ is (Ker Li) ± , the L 2 orthogonal 
complement of the linear span of the {9jipo}j—i- 

Proof. We proceed similarly to |Kik08l Lemma 5.4.3]. The n = ground state solution tpo to ( 13.4b is characterized 
in IILie77l as the solution, unique up to translation and phase rotation, to the following minimization problem: 

E(<po) = I„ ■= inf{S(0); e H 1 (M 3 ) 1 ||<^|| 2 2 = p}, (3.8) 

for certain p > 0. We claim that this implies that L\ > on {(po} ± . Indeed, let \\v \\^2 = ||<^ ||l 2 ; (v, fo) — 0- F° r 
s G ( — 1, 1), define <f) s = (1 — s 2 ) 1 / 2 ip + sv, so that Q(cf> s ) = Q(fo)- Calculating that (f> s \ s=0 = Vo. 9 S \ S=0 4> S = v, 
d 2 \ B=0 (f> s = —(fo we deduce from ( 13.81 ), 

< d 2 \ s=0 E(cf> s ) = (E'(Lp ),-Lp ) + (E"(tp )v,v) = -LV (Q'(Lp ),(po) + (E"(cp )v,v) = (v, (E" -oj Q")v), 



establishing the claim. We took into account that ipo satisfies the stationary equation E'(cpo) = looQ'(<Po) an d a ls° 

\'h = nv\\h 



that (Q'(<p ),<p ) = 2\\Lp \\ 2 = 2\\v\\ 2 = (Q'(v),v) = (Q"v, v) 
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So L\ is non-negative on a codimension one subspace. On the other hand, since the integral kernel of A -1 is 
strictly negative, while ipo is strictly positive and Lo(po = 0, it follows that ((poL\ (fo) < so that there certainly 
exists one negative eigenvalue characterized as 

-A := M{{v,L lV ); \\v\\ L 2 = 1} < 0. 

Let rjo be the corresponding eigenfunction, Lir/o = — Ao??o- To prove that (— Ao,0) C p{L\), the resolvent set, 
consider the minimization problem 

inf-C^iiw); ||«|| a = 1, (rn,v)=0}. (3.9) 

Now the relation Lqlpq = 0, together with translation invariance, implies that Lidjipo = 0. Moreover, it is proved in 
IILen09l that (fo is nondegenerate, in the sense that the kernel of L\ is spanned by the djipo, 1 < j < 3. Hence, by 
consideration of linear combinations of the eigenfunctions 770 and dj ipo, that the number defined by d3.9t is < 0. In fact 
it must equal zero since if it were negative a simple compactness argument (based on the negativity of Uq) would imply 
the existence of a negative eigenvalue in the interval (— Ao, 0) with corresponding eigenfunction 771 orthogonal to 770. 
But since 770 , 771 would then be an orthogonal pair of eigenfunctions of L\ with negative eigenvalues, and both having 
non-zero inner product with i^q, this would immediately contradict the fact that L\ is non-negative on {(po}^- □ 



We conclude with a few remarks on the stability of solitary waves to the Choquard equation. By Remark I3T31 we 
know the w-dependence of a localized solution 4> LJ (x)e~' tu ' t to d3.3t : one has <j) u {x) — ( 2 u((\x\), where £ = \J-lo. 
From this we can obtain the frequency dependence of the charge: 

Q(u) = f <f> 2 (x)dx = ( i j u 2 {(x)dx = ( [ u 2 {y)dy = {-uj) 1 ' 2 [ u 2 {y)dy. 



It follows that for all negative frequencies ^ < . By the Vakhitov-Kolokolov stability criterion ( IIVK73I ). this leads 
us to expect the linear stability of no-node solitary waves (the ground states) in the Choquard equation. 

Proposition 3.5. The ground state solitary wave tpo(x)e~' lUot of the Choquard equation ( 13.31 ) is linearly stable. 



To determine the point spectrum of 



L a 
-L 1 



observe that if 



is an eigenfunction corresponding to the 



eigenvalue A G C, then —X 2 R = LqL\R. If A 7^ 0, then one concludes that R is orthogonal to Ker Lq = {(fo}, hence 
we can apply Lq 1 ; taking then the inner product with R, we deduce that: 

-X 2 (R,L Q 1 R) = (R, LxR) , 

which implies that A 2 6 K. Moreover, by (13 .91 . A 2 < 0, leading to the conclusion that the point spectrum ad(JL) C 
iR, and hence the absence of growing modes at the linearized level. The (nonlinear) orbital stability of the ground 
state solitary wave was proved in IICL82I . 

Remark 3.6. In view of 1CGG121IBC12I . one expects that the linear stability or instability of small amplitude solitary 
waves is directly related to the linear stability or instability of the corresponding nonrelativistic limit, which for Dirac- 
Maxwell is given by the Choquard equation. We hope that this may provide a route to understanding stability of small 
solitary waves solutions for the Dirac-Maxwell system. 



4 Proof of existence of solitary waves in Dirac-Maxwell system 



In this section, we complete the proof of Theorem 11.11 It is obtained as a consequence of Proposition 14.61 after the 
application of a rescaling motivated by the discussion in $3] 



We write d>(x) 



4>i{x) 
(j)2{x) 



, where for j = 1, 2 the <fij e C 2 are essentially the components of <f> in the range of the 



projection operators LIi = ^(1 + /?), and 77 2 = i(l — j3) (under obvious isomorphisms of these subspaces with C 2 ). 
Applying LI\ and LI2 to (11.61) . we have: 

u<f>i = <r-(-iV - eA)(j> 2 + mcpx + eA ^, (4.1) 

tufa = a-{-iV - eA)0i - to02 + eA> 2 , (4.2) 
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- AA° = e{4>l4>i + ^0 2 ), - AA = e4*(x<t> = e(^<r0 2 + fatrfa) . (4.3) 

We write (14.3b as 

A° = eN* (<ft<h + <&<h) , A = eN*(^ 2 + ^i), (4.4) 



and regard the potentials A and A = (A 7 ) as non-local functionals of 1 



Newtonian potential. In abstract terms, the equations are of the form u)Q! = £' where the charge functional is 



Above, N(sg) = (47r|x|) _1 is the 



Q0) = / <j>*{x)<j){x)dx, (4.5) 
and, regarding A , A as fixed non-local functionals j4Aj of tfi, the Hamiltonian £(<fi) is given by 



For future reference we recall the following trick from |Stu99|: 

Lemma 4.1. Let £ Q be a finite collection of vector fields on the phase space which are infinitesimal symmetries, in 
the sense that (Q' , £ Q ) = = (£' , £ Q ). Then any solution of the equation ujQ' — £' — a a t; a = ,/or some set 
a a G R, is also a solution ofcuQ' — £' = 0, as Zong as f/ie matrix £w is we/Z defined and nondegenerate. 

Proof. For sufficiently regular £^ it is possible to take the inner product, yielding ^2 o, a C' 3 ) = which gives the 
result. (The precise meaning of sufficiently regular is just that this computation is valid; it would be sufficient for £ Q 
to he in a subspace F of L 2 with the property that the equation luQ' — £' — ^ a a ^° = holds in the dual of F.) □ 

Example 4.2. For tp : K — > C and Q = | J |?/>| 2 and £ = J ^\Wtp\ 2 — ^jjrj-|'0| p+1 the symmetry of phase rotation 
corresponds to the infinitesimal symmetry = iip, and it is easy to check that given an H 1 distributional solution 
of ojQ' — £' — a£ = 0, i.e. a weak solution of — — \ip\ p ip = cutp — iaip, for any a G K, one necessarily has a = 0. 
The same holds in higher dimensions as long as p is such that the equation holds as an equality in H . 

Remark 4.3. The advantage of solving the more general equation with the unknown "multipliers" a a is that in an 
implicit function theorem setting the multipliers can be varied to fill out the part of the cokernel corresponding to the 
symmetries. It is then shown after the fact that the multipliers are in fact zero. The choice of is determined by 
the symmetry group; in the case of Dirac-Maxwell the relevant group is the seven dimensional group generated by 
translations, rotations and phase rotation. The infinitesimal versions of these actions give the following vector fields 
( IIBD641 ): 



£ = V0, T)=- 



er 

CT 



C = i<j>- (4.7) 



In accordance with the heuristics in $2] we introduce functions •P 1 (y,e), ^2(2/: e ) £ C 2 and A^(y,e) by the 
following scaling relations: 

4>i(x,lj) = e 3 <£i(ex,e), cf> 2 (x, u) = e 2 $ 2 {ex, e), eA°(x, u) = e 2 A°(ex, e), eA j (x, w) = e 3 A j (ex, e), (4.8) 



where e and ui are related by u = — V m 2 — e 2 . Then, writing Vj, for the gradient with respect to y 3 = ex 1 , 1 < j < 3, 
we have: 

- 2m$i + ier-V y <P 2 - e 2 A°<P 1 = -(m + cj)#i - e 2 A<r<£ 2 , (4.9) 
1 -<P 2 + ia--Vy<P 1 - A°<P 2 = (— — )&2 - e 2 A-cr<P 1 , (4. 10) 



2m ' y 1 ' \2m 

A = e 2 N* ($* 2 $ 2 + e 2 $* 1 $ 1 ^) A = e 2 N * ($la$ 2 + <P 2 er$^) . (4.11) 
Let ipo G ^(M 3 ) be the ground state solution to the Choquard equation with ujq = — t^-: 

-^^o-7^A^ -(N*^ 2 )^o = 0. (4.12) 
2m 2m 
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That is, ipo is a strictly positive, spherically symmetric, smooth, and exponential decaying function. As discussed in 
the previous section, such a solution exists by IILie77l ; the value ujq = — (2m) _1 is chosen for our convenience. Using 
ipo, we can produce a solution to (l4.9b -( t4~TTb in the nonrelativistic limit e = 0: 



& 2 



with <E>2 = 



$i = —(T-Vy<P 2 : 

2m 



A° = e 2 N*^ 2 , A 1 = E*(p Q d 2 <po, A 2 = +— N * ip diip , 

m m 

The symmetry of this configuration is axial, with the magnetic field along the z axis of symmetry. 



A d = 0. 



(4.13) 
(4.14) 



Lemma 4.4. Let<P 



e ii^R 3 , C 4 ). Then A^ defined by gTTJ satisfy 



A' 1 e L°°(R ), < n < 3. 



Proof. The functions A M defined by (14.11b are of the form N * h with h := fg, where f,g£ i? 1 (R 3 ). Due to 
the Sobolev embedding iJ^R 3 ) C L 6 (R 3 ), we have h £ L P (R 3 ) , 1 < p < 3. By the Holder inequality, one has 
\(M~ 1 Xb 1 (x)) * h\ < \\\x\~ 1 \\ L 3/2(B 1 )\\ h \\L s < oo, |(1 - XBi(x)) * h\ < \\h\\ L i < oo, where B 1 is the unit ball 
in R 3 and \b 1 is its characteristic function, hence * h € L°°(R 3 ). Furthermore the structure of ( 14.1 It makes it 
clear that the mappings e) h-> A M , < /i < 3, are smooth from iJ^R 3 , C 4 ) x R to L°°(R 3 ). □ 



Define 



X = H\M. 3 ; C 2 ) © i? 2 (R 3 ; C 2 ) C L 2 (R 3 ; C 2 ) © L 2 (R 3 ; C 2 ) 
Y = H^M 3 ; C 2 ) © i 2 (R 3 ; C 2 ) C L 2 (R 3 ; C 2 ) © L 2 (R 3 ; C 2 ). 



(4.15) 
(4.16) 



Introducing f= — ier-'Vy and substituting u) = —y/m 2 — e 2 , we rewrite ( 14.9b . (14.101 as the equation T = 0, where 



J 7 : X x (—m, +m) 



r. 



J" 



#1 

# 2 



2m#i + + e 2 A <?i - (to - \fm 2 - e 2 )<£i - e 2 A-<x<2> 2 
-L# 2 + ^ 1+ A # 2 + fi 



1 



2m 



\2m 



+ yjm 2 - e 2 



<5o 



£ A-CT^! 



(4.17) 



As above, we regard the A^ = (A , A), A = (A- 7 ), as non-local functionals A M = A M (<£, e) determined by ( 14.111 ). 



With this understood, the entire system is encapsulated in the equation J-^, e) = for $ = 
the original variables: 



<J>o 











only. In terms of 



(4.18) 



where the functionals Q, £ are defined by (14.51 . (14.6b . The nonrelativistic limit satisfies J"(<Z>, 0) = (cf. ( 14.13b . ( 14.14b ). 
so that to obtain solutions for small e it is necessary to compute the derivative of T at the point (<P, 0). This is de- 

"ll 



termined by the set of directional derivatives. Let e\ 



and e 2 



, and let g g H 1 ( 



To compute 



the directional derivatives first note that A 3 drops out on putting e = 0, and then note further that by ( 14.1 lb only the 
derivative of A at (<?i , 0) with respect to <S> 2 is nonzero, with derivative given by 



dt 



A" 



( U + J ' e )l*=o.-o = 2e 2 N* (^ Re(e 1)5 ) c2 ). 



We deduce that for C 2 -valued functions U and V, 



dt 



I 



&2+tV 



, c 



M 



where 



M 



2m f 
f - 2 ^+A + 2ip e 1 N*(i Po Be(e 1 , 



-icr- V„, 



N 



47r|y| 



(4.19) 



Thus, the derivative of T at the nonrelativistic limit point (<P, 0) is the linear map DJ-($, 0) given by the matrix M. 



This is a differential operator, which we consider as an unbounded operator on L (I 



L 2 



5 ;C 2 ). 



Lemma 4.5. 



'u 




F 


V 




G 



7. 77ie map M 

2. M fflapi X continuously into Y (cf. ( 14.15b , d4.16l >). 

3. 77ze kernel of M is g/ven fey 



is a Hennitian operator with domain X. 



KcrM = 

4. The range of M 

Range M = (KerM)" 



{ ( ~ 2m ' V ) ' V = a ' V » Vo ei + ib<f0 ei + c( ^o e 2 , a G E 3 , 6 G E, c G c| 



"[/" 




F 


V 


1 — ^ 


G 



is closed in the topology ofY and is given by 





F 


{ 


G 



£Y:Be('—-G)£(KerL 1 y 



Im 



(IL 

\2m 



G G(KerLo)- 1 -, 



V 2m 



G G (Kcr L )- 



where means the orthogonal complement with respect to the inner product in L 2 L 2 



5. The inverse of M 



is given by 



2m v 



V = e 1 Vi + e 2 V 2 = Lr 1 R° 



V 2m 



G 



iL 1 Im 



(tL- 

\ 2m 



G 



where the definitions and properties of the operators Lq, Li are given in 



Proof. The proof depends on some properties of the linearized Choquard equation from [Len09| which are stated in 
Sj3] The fact in (1) that M is Hermitian follows from the fact that pis Hermitian. From Lemma [4~4l the assertion (2) 
is immediate from the properties of N and the fact that <po and its partial derivatives are smooth and exponentially 



decreasing. To prove (3), (4) and (5) we consider how to solve M 



r 




F 


V 




G 



, i.e. the system 



M 







V 





2mU + fV 
A°V + 2^ eiN* (ipoBeVi] 





F 




G 



We first express U in terms of V by U = ^(F — fV) , and, writing V = Vid + Vi^-i. 

ip jp A V V 

jr- + - — +A°V + 2u oei N * (w Re Vi) = G. 

Zm im im 

Referring to the definitions in $3]of Lq and L\, with ojq set equal to — (2m)" 1 , we arrive at the following equations: 



_ f fF_ 



\2r 



G 



L V 2 



\2m 



G 



(4.20) 



It is useful here that the components with respect to ei and e 2 are decoupled. Noting also from the form of Lq, Li that 
these operators take real/imaginary valued functions to real/imaginary valued functions, and further that L\ = Lq on 
pure imaginary functions, we obtain the given formula for V, and hence for U, immediately from Sj3] The identification 
of the kernel in (3) is then a specialization of this, given the information on Kcr Lo and Kcr L\ in SJ3J and also (4) is a 
consequence of the identification of the ranges of Lq and L\ given in <J3] □ 
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The statement of Theorem ll.il will follow from the following result. 

Proposition 4.6. There is e* > such that for e G (— e», e*) there is a solution to ( 14.1l l-( l4~3"1 ), w/f/z w 
g;ven by the ansatz d4.8t w/f/z (<?Xi ^2) obtained as the image of a C°° -function 



-ym 2 



$er((-£»,e«); (^(R'jC 2 )© i? 2 (M 3 ;C 2 )) nKerM) 

satisfying <Z>(0) = 4 and wirt A G C*°° ((-e*, e*) ; 1 n A-? G C*°° ((-e*, e*) ; i? 1 n given fry (HTTTT i: 

A = e 2 N * + 6 2 ^<£ 2 ) G C°° ((-e„ e.) ; tf 1 n , 

A = e 2 N * (<^CT# 2 + *5o^i) G C°°((-e*, e*) ; i? 1 n . 

Above, H 1 = i? 1 (IR 3 ,K) is the homogeneous Dirichlet space of L 6 functions with \\f\\ 2 ^ 1 = / |V/| 2 dx < 00. For 
smaW e, one has 

(4.21) 



|*2-*a||ffi =0(e 2 ) 



77ie are even in e. 



Proof. Solutions of d4.U - d43b for small e can be produced by solving T = 0. The proof of existence of solutions 
to this equation is by the implicit function theorem and Lemma |4~T1 perturbing from the nonrelativistic limit point 
J-($, 0) = 0. To start we claim that F, as defined in ( 14.171 ), is a C°° function X x (— m, +rn) —> Y. To prove 
this notice that the expression for T is manifestly smooth in e for e 2 < m 2 , and its dependence on <Pj is built 
up from compositions of certain multilinear maps and linear operators; the structure of the expressions obtained after 
successive differentiation is the same. Referring to the specific formulae, the fact that these expressions are all C°° is an 
immediate consequence of the fact that multiplication gives continuous bilinear ( smooth) maps H 1 x H 2 —> H 1 
and H 2 x H 2 — > H 2 (Moser inequalities) and Lemma [4~4l We are looking for <P(e) in the form 



<P(e) =$ + &(e), <F(0) = 0. 



(4.22) 



We use the same component notation as above: 3 
the function 

G: Xx 

g(^,a,b,e) = + e) 



We apply the implicit function theorem to 



xl 3 x ( — m, +m) 



a-V, ; 



Y, 



"e(#i +!?!)" 




ecr(£>i - 


f 






a(& 2 H 


-#2). 



(4.23) 



Remark 4.7. Referring to ( 14.31 ) we have introduced a linear combination of the six infinitesimal symmetries corre- 
sponding to translation and rotation. The action of phase rotation is not independent of rotation in the nonrelativistic 
limit, which is why the seventh parameter does not appear. In terms of the original variables (cf. ( 14.71 )): 



g(^a,b,e) = 







£' - ujQ! + ea • £ + e 2 b • 77 ) evaluated at </> = 



e 2 (&2 



#2). 



Computing the derivatives of ( 14.23b at e = 0, we see that the linear span ({d a Q,db G', 1 < .7 < 3}) is equal to 
Ker M. Referring to Lemma |431 this establishes that the derivative of Q at e = 0, \P = 0, a = 0, b = with respect 
to ('J, a, b) is a linear homeomorphism from ((Ker M) 1 - n X) x I 3 x I 3 onto Y. It follows that there is e* > such 
that there exist C°° functions e i-» (^(e), a(e), b(e)) elxl 3 



g(^(e),a(e) ) b(e),e) =0, 



R 3 , defined for e G (— e*, e*), such that 
<p-(e) _L KerM, e G (-e*,e»). 



(4.24) 



This latter condition serves to divide out by the action of the symmetry group, giving a local slice. 

Referring to Lemma |4~T1 to deduce that these in fact generate solutions of T = 0, for sufficiently small e > 0, it 
is sufficient to verify that a(e) = 0, b(e) = 0, which is in turn a consequence of the nondegeneracy of the matrix of 
inner products of the infinitesimal vector fields, scaled as above. This amounts to the need to verify nondegeneracy of 
the 6x6 matrix 

~{d y j(t>,d y k4>) (d yj (j),%'E k >(f>) 



for small e. (In the matrix ( 14.25b the indices j, j', k, k' run between 1 and 3.) 



(4.25) 
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Lemma 4.8. The matrix given by ( I4.25I ). evaluated at • 







V 2 (S 2 - 





is nondegene rate for small e. 



y=ex 



Proof. Clearly the dominant terms arise from the second ("large") component giving rise to diagonal matrix elements 
which, referring to the block form in ( 14.251 ). are 0(e 4 ). Since \Pj = 0(e), the result will follow from nondegeneracy 



of the matrix with <Jj set equal to zero. Using e 
term: 

A 



$2 



and <S>-} 







(d y j(j),d y k(j)) 



dyjipo, (- 



4m 2 



, we calculate the first diagonal 



)d y k<p Q ) + (d y j<po, d y k<po) + 0(e) = -^-(^ipo,(-A y )ipoj + 0(e), 



where we took into account the spherical symmetry of ipo, which leads to (d y iipo,dyiipo) = ^(fo, (~^-y) l Po)- Next 
for the off-diagonal terms we compute, again using the same expression for e~ 2 0: 



e~ A (d y j<t>, -E fc ^) = 



4m z 



dyjf 



tpo 




d yJ 



ipa 




2 CTfe ' 



ipo 




) +0(e) = 0(e), 



The first two terms are identically zero since ipo is spherically symmetric (so that by parity considerations it is L 2 
orthogonal to all of its first partial derivatives, which are in turn orthogonal to all of the second partial derivatives). 
Finally, for the second diagonal term: 



6j>k' 



(ipo,ipo) +0(e). 



The non-degeneracy of the matrix (14.25b for small e follows. □ 
Since the implicit function theorem proves that <Pi(y, e), ^2(2/, e) are C°°-functions of e G (— e*, e*), one has 

||#i-# 1 || ir »=0(e) ) \\$ 2 - £ 2 \\ H i = 0(e) . (4.26) 
To prove a stronger estimate (14. 2U . we take the derivative of ( 14.17b with respect to e at e = 0; this yields 

Md^\ c=a = 0, 

with M given by (14.19b . Due to ( 14.22b , one has d e <P\ e=a = d e &\ e=0 ; and the requirement ( 14.241 ) leads to d e <P\ e=0 = 
9 e !f r | €=0 = 0, and hence 

||<2>! - ^11^2 = 0(e 2 ), \\<P 2 - S 2 \\ Hi = 0(e 2 ). 
Finally, notice that since the explicit dependence of T is on e 2 , we have 



F(^(-e),^ 2 (-e),-e) = ^(^(-e),^ 2 (-e),+e) =0 



(4.27) 



and hence, <Pj (e) = <Pj (— e) since otherwise it would be possible to contradict the local uniqueness part of the conclu- 
sion of the implicit function theorem (applied to Q with a = 0, b = 0). This completes the proof of Proposition 14.61 
and thus of Theorem ll.il □ 

Remark 4.9. The solutions of T = are obtained for both positive and negative epsilon close to zero, but the e negative 
branch apparently gives rise to solutions of the Maxwell-Dirac system via d4.8b which are related to the positive branch 
as follows. By ( 14.8b , the branch which corresponds to negative e has the form 



4>(x,w) 



<k(x,w) 




(j> 2 (x,U))_ 





(-e) 3 0i(-eir,-e)' 
(-e) 2 $ 2 (~ex,~e) 



-e 3 <P 1 (-ex, e) 
e 2 $ 2 (-ex,e) 



-\J m 2 



e 2 , 



e > 0, 



where we took into account that e), < P 2 (y, e) obtained from Proposition 14.61 are even in e. Comparing to ( |4.8b , 
we conclude that this branch is related to the e-positive branch <fi(x, w) by 

<f>i(x,oj) = -4>i(-x,u), 4>2(x,w) = (j> 2 (-x,u), 

so that A°(x) = A a (— x), A(x) = —A(—x). That is, these two branches have the same magnetic field but opposite 
electric field (see 1BD64I §2.3 and §5.4]). 
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Remark 4.10. We briefly consider the symmetry properties of the solitary wave solutions: in HLis95l Section 5] Lisi 
gives an ansatz for the solitary waves, using cylindrical coordinates (p, z,8), from which symmetry properties can be 
deduced. For our situation the relevant ansatz for the Dirac wave function is 



(4.28) 



It seems likely that the solutions constructed via Proposition l4.6l have this symmetry and that this fact could be proved 
via an application of the implicit function theorem within the symmetry class of ( 14.281 ). 
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